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ESTIMATES FOR CERTAIN INTEGRALS OF PRODUCTS OF SIX 

BESSEL FUNCTIONS. 

DIOGO OLIVEIRA E SILVA AND CHRISTOPH THIELE 


Abstract. We establish good numerical estimates for a certain class of integrals involv¬ 
ing sixfold products of Bessel functions. We use relatively elementary methods. The 
estimates will be used in the study of a sharp Fourier restriction inequality on the circle 

in p. 


1. Introduction 

Let (S^, cr) denote the unit circle in the plane equipped with its arc length measure. The 
companion paper [T] discusses partial progress towards understanding the optimal constant 
Copt in the endpoint Tomas-Stein adjoint restriction inequality [6] on the circle: 

(1) II/<^IIl6(m 2) < Copt II/IIl2(§i), 

where the Fourier transform of the measure fa is given by 

(2) fcr{x)= [ f{uj)e~'-^'‘^daui, (x G M^). 

Jsi 

It is conjectured that equality is attained in 0 when / is a constant function. For the 
constant function / = 1 , the sixth power of the left-hand side of inequality 0 turns into 
the integral 


poo 

(3) (27r)^ / Jo{r)rdr, 

Jo 

where the Bessel function of order n, denoted Jn, is defined via the identity 

(4) = 27r(-i)”J„(|x|)e*’^”§(*). 

Part of the analysis in [T] consists of a Fourier expansion of / on the circle, and one needs 
estimates with rather precise numerical error bounds for integrals of the following two 
types: 
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( 5 ) 

and 


In =/( 


0,m,n 



Jn-\-m {r)Jn{r)Jm{r)Jl{r)rdr 


POO 

(6) h = h,m,n := / Jn+ni{r)Jn{r)Jm.{r)Jl{r)Jo{r)rdr. 

Jo 

The purpose of the present paper is to establish these estimates, summarized in the fol¬ 
lowing theorem: 


Theorem 1. Let n >2 be some integer. Then eaeh of the following quantities is less than 
0.002n-^; 


(i) For n > 7: 


and for n >3: 


(ii) For any n > 2 : 


^ 3 13 

-'0,0,n--1- 


1 


n 32 'k‘^ {n — l)n{n -t- 1) 

1 


T _J_i_ 3 

^1.0,n .0 „ 


dvr^ n 3271^ (n — l)n(n -I- 1) 

15_1 

2 n{n -|- 1 

1 


Gdvr^ n(n-|-l)(n-|-2) 


h,2,n — 


9 


GdTT^ n{n -|- l)(n -|- 2) 
Moreover, each of the following quantities is less than O.OOlSn”"^: 
{in) For n > d; 

^ 1557 1 

o,4,n X02d7r2 n(n -|- 1) (n -|- 2) (n -|- 3) (n -|- d) 


7l,4,n 


855 


1 


102d7r2 n(n -|- l)(n -|- 2)(n -|- 3)(n -|- d) 
{iv) For even m > 6 and n > m: |/o,m,n| o.nd 


Thus Theorem controls integrals of the two types Iq and Ji for n > 2 and even 
0 < m < n, with the exception of the five cases m = 0 and n = 2,3, d, 5, 6 for Iq, and the 
two cases m = 0 and n = 2,3 for Ii. It follows from Tablebelow that, in these exceptional 
cases, the quantities are still less than O.Oln”"^, which provides information about Iq and 
Ii with at least two percent relative accuracy. 

It follows from the methods of this paper, or alternatively from general principles, that 
such a result holds with bounds cn“^ for any positive number c in place of 0.002 or 0.0015, 
and with some finite set of exceptions. The point of Theorem is to narrow down these 
exceptions precisely for the specific numbers c = 0.002 (for m = 0,2) and c = 0.0015 (for 
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m > 4). Slightly better numerical estimates are listed in Sections and for the various 
cases, but for simplicity we do not reproduce all of them here. 

Our methods apply to obtain a more general set of estimates than the ones listed in 
Theorem but we focus on the stated estimates which are needed in [T] . There exists 
a very satisfactory theory of similar integrals of products of two Bessel functions, see for 
example Lemmata and below, and a still explicit but substantially more complicated 
theory for integrals of products of four Bessel functions. While integrals of sixfold products 
of Bessel functions still fall into the class of functions for which explicit symbols have been 
introduced in the theory of hypergeometric functions and their generalizations, we do not 
know how to obtain our rather accurate numerical bounds in a much easier way than by 
the elementary but somewhat laborious approach presented in this paper. 

Our approach is to expand four of the six Bessel function factors, namely those four 
with the lowest orders, into their asymptotic expansions. This will reduce the integrals in 
question to core integrals of the type 


roc roc 

(7) / Jnir)Jn+m{r)sm{ir)r~’‘dr, / Jn{r)Jn+m{r)cos{(.r)r~^dr 

Jo Jo 

for i = 0, ±2, ±4. For these integrals, one has good information as in Lemmata EE 
In more detail, the paper is organized as follows. In Section we review the theory 
of Bessel functions inasmuch as it is useful for our purposes. In particular, we establish 
the aforementioned lemmata, together with asymptotic expansions with precise control on 
the error terms. In Section we prove some useful estimates for binomial coefficients, 
the Gamma function, and the coefficients that arise in the various asymptotic expansions. 
The analytic part of the proof of Theorem begins in Section where we asymptotically 
expand the functions Jq and Ji. Section accomplishes the same for the function of next 
lowest order, namely Jm- Finally, Sectionis devoted to the analysis of the core integrals. 
The estimates from Sections EHH are then assembled together in Section The approach 
works for n > 20, and so for n < 20 we numerically estimate the integrals; this is the 
content of the hnal Section [8l 

We close this discussion with a brief illustration of the difficulty involved. Figure 
depicts the plot of the integrand of /i, 6,9 between r = 0 and r = 100. One observes an 
initial region until about r = n = 9 where the function is very small. Then one sees a 
region with fairly erratic behaviour until about r = n? = 81. Past r = 81, one sees a more 
repetitive behaviour where one has good asymptotic control. The asymptotic region yields 
a positive contribution to the desired integral, which is in general of the order n“^. The 
erratic region yields a negative contribution which nearly cancels the positive part from 
the asymptotic region. In question is a very good numerical control of the order n~^ of the 
small difference. The main tools to capture this cancellation are the algebraic identities 
from Lemma and an exact orthogonality formula due to Kapteyn [2] which can be found 
in Lemma [3l 
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2. Background on Bessel functions 

We rewrite the definition Q of the Bessel function in the form of a Bessel integral which 
is the starting point in 0 p. 338]. For n G Z and z > 0, we claim that 

( 8 ) Jn{z) = — r 

J-n 

More precisely, replacing 0 = w -|- 7r/2 in Q and using even symmetry of the cosine we 
obtain for the right-hand side of Q: 

(_j)n fir n . 

(9) cos(na;)(ia;, 

2^ J-n ^ Jo 

from which the equivalence of Q and Q is evident. 
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The Bessel function, defined via Q for general z G C, is an entire function. From Q 
we obtain the estimate 

(10) |Jn(^)| 

Differentiation under the integral sign in Q and integration by parts yield the following 
recurrence relations in the sense of meromorphic functions: 

(11) Jn-l{z) - Jn+l{z) = 2j'^{z), 

2/7- 

(12) Jn-l{z) + Jn+l{z) = —Jn{z). 

A different representation of the Bessel function is the Poisson integral, which contains 
a power of a trigonometric function rather than a power of an exponential function: 

(13) Jn{z) = ^/2)r(l/2) Jo cos(zcos(0))sin2^(0)d0 . 

Here we use the Gamma function 

roo 

r(s) := / , 

Jo 

which satisfies the functional equation sr(s) = r(s + 1) and thus meromorphically extends 
the factorial, that is r(n + 1) = n! for natural numbers n. We mainly need the Gamma 
function for half-integer values, which can be expressed as 



This can be recursively verihed from r( 2 ) = which in turn can be read from the 
well-known property 


(15) 


sin(7rx)r(x)r(l — x) = vr. 


The latter can be seen by verifying periodicity of the left-hand side together with growth 
estimates which force the left-hand side to be constant. 

, one verifies the 


To see equivalence of the Poisson integral representation (13) with 


case n = 0 by substitution and then verifies by partial integration that the Poisson integral 
also satisfies the recursion relations © and ( |12[ ). Combining these two second order 
recurrence relations into a first order relation between and Jn+i, equivalence of the two 
integral representations follows recursively by a uniqueness result for ordinary differential 
equations, where we use that both integral representations vanish for z = 0 and n > 0. 

From the Poisson integral representation one sees the following estimate from [3 §3.31(1), 
p. 49], useful for small 2 :: 


\J'n{z)\ < 


|^|ng|S2| 

2’^n! 


(16) 
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where we have used 


sin^”(0) cie 


n! r(n + l/2)r(l/2) 7o 

which one proves by induction on n using integration by parts. 

We turn to the core integrals Q. The case £ = ±2 will be the most pleasant to deal 
with via the following lemma: 

Lemma 2. Let 0 < n,m have the same parity, and 1 < k < n + m. Then if k is even 

poo 

/ Jn{r)Jm,{r)r~^ cos(2r) dr = 0 , 

Jo 

and if k is odd 

roo 

Jn{i')Jmir)r~^ sin{2r) dr = 0 . 


f 


f 


Proof. By the parity assumption, we may extend the integrals to the full real line. It then 
suffices to show that the Fourier transform 

r*00 

Jn{r)Jm{r)r~^e~‘^^'^ dr 

' —OO 

vanishes at ^ = 2. Substituting ^ = cos a; on the right-hand side of ([^ yields 

7 _ (-*)” r j^£cos(narccos0 

Mz)- ^ J_^e df. 

Hence we see that Jn is the Fourier transform of the function 

BniO = , 

TT i ’ J 

where T„ denotes the Chebysheff polynomial T„(^) := cos(n arccos 

We first consider the case k < n in the lemma. Then the Fourier transform B 
Jn{f)r~^ is still supported on [—1,1] since Bn has vanishing moments of orders 0,1,..., A: — 
1. This can be deduced from (13). Seen as the convolution of an IP function with an IP 
function for p = 2 , the function Bn ^ * Bm is continuous. Since it is also supported on 
the interval [—2,2], it must vanish at ^ = 2. This proves the lemma in case k < n. If 
n < k < n + m, we distribute some powers of r over Jm and argue similarly. □ 

The understanding of the dominant case £ = 0 of the core integrals ([^ begins with 
Kapteyn’s identity, proved in a delightful two page paper [2]. 

Lemma 3 ([2]). If n,m > 0 and n + m ^ 0, then 

poo 

( 17 ) / Jn{r)Jmir)r~^dr = - 


(-k) r 

n of 


2 sin^^TT 


TT 


with the following natural interpretation in case n = m: 

1 


f 


Jniryr ^dr = 


2n 
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Note in particular that © vanishes if m — n is a nonzero even integer. Moreover, 
identities (14), (15), and some algebra yield 


2 sin ■ 


-TT 


2-ir(™±n) 


TT — ri^ 


n 




)r( 


n—m-\-2 


)r( 


m—n+2 ' 


2 \ 2 \ 2 / 
From Kapteyn’s identity, one obtains the following more general result. 


Lemma 4. If 0 < n,m and 1 < k < n 

roo 

(18) / Jn{r)Jm{r)r~^dr = 


m, then 

2-^r(A:)r( 


m+n+1—> 
2 > 


r(i 


m+n+l+/c n—m+fc+1 \ f m—n^k-\-l 


)r(^ 


2 V 2 y 2 

Proof. The identity is true for A: = 1 by Lemma Let k > 1 he given and assume that 
identity ( [I^ is true for this particular k. To prove the identity for /c + 1 we denote the 
integral in (|l8[) by Tn,m,k- Then, by the second recursion in (12) and using the induction 
hypothesis, we have that 

^n,m,k+l — ~h ^n+l,m,k') 


2-(fc+i)^-ir(fc)r( "^+^-^ ) 


p^ m+n+fe ^p^ 7i—m+fc ^p^ 

2-(fc+i)r(A: + l)r( ""+^-^ ) 


m—n+k-\-2 

2 


+ 


r( 


2-(fc+i)^-ir(fc)r( "^+^+^-^ ) 

m—n+k > 
2 > 


m+n+2+/c ^p^ n—mH-fcH-2 


)r( 


2 y 2 
(m + n + k){n — m + k) + {m + n — k) {m — n-\-k) 


r( 


m-\-n-\-2-\-k 


)r( 


n—m-\-k-\-2 


)r( 


m—n-\-k+2 


‘ink 


2 V 2 y 2 ) 

The second fraction in the last line is equal to 1, and this proves the inductive step. □ 


Note in particular that if k is odd and satisfies k < |n —m|, then (18) vanishes since the 


function l/T has zeros at s = 0, — 1, —2,... Note also that three of the Gamma factors, the 
one with argument k in the numerator and the two involving the difference m — n in the 
denominator, typically form a binomial coefficient, which can be roughly estimated by 2^. 
An alternative approach to Lemma is the integration theory of Weber and Schafheitlin 
as outlined in [TJ §13.24, p. 398]. 

The case (. = ±4 of the core integrals ([^ gives small error terms which we estimate with 
the following lemma. 


Lemma 5. Let 0 < n,m and 1 < k < n + m. Then: 


f 


Jnir)Jm{r)r-’^e^^^dr 


')k—l 


< 


(n + m — A:)! 


^n+m 


Proof. We estimate this integral by the descent method, changing the contour of integration 
to the contour consisting of a line segment from 0 to iN for some large N, then a line 
segment from iN to N, and then a ray from to oo along the real axis. Only the first 
integral provides a substantial contribution, since the next two segments give contributions 
that tend to 0 as tends to infinity. 


Along the first segment of the contour the integral can be estimated using (16) by 
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^OO -j 1 /•OO 

2^-1 _|_ 

“ 4n+m 


The integral over the second part of the contour is estimated using (10) by 


j-N 

\/2 / \Jnix + i{N - x))\\Jm,ix + i{N 

Jo 


x))\\x + i{N - x)\-^e-<^-^Ux 


rN 

<V2 {N/V2)-^e-^^^-^'>dx < 2^+^Ar-^ , 

Jo 

which tends to 0 as —>■ oo. The integral over the third piece of the contour is estimated 
using the fact that the functions Jn and Jm are in as Fourier transforms of func¬ 
tions and thus the continuous function is in L^. It then follows from the 

dominated convergence theorem that 


roo 

lim / \Jn{r)Jn+m{r)\r~^dr = 0 . 

N^oo Jj^ 

Adding the contour integrals and letting N ^ oo yields the desired bound. □ 

To arrive at the core integrals, we need asymptotic expansions of the Bessel functions 
near infinity as in the following lemma. 

Lemma 6. Let n G N and 51?(2;) > 0. Let ujn = z — n-nI2 — tt/A. Let i € N be sueh that 
i > max{n — 1/2,1}. If i is even, then 


(19) 


- —1 
2 


'Jn{z)= [ — ) (cOSWn) ( y~l(-l) 

. k=0 


k O-2k[n) 

■y2k 


k=0 




If i is odd, then 

( 20 ) 


2 


'£-3 


(cOSWn) j y~](-l) 


, k=0 


k <^2k{n ) 

z2k 


(sinwn) 


ka2k+i{n) 


,k=0 


z 


2k+l 


+ Rn,l{z) 


Here the eoefficients aj{n) are defined by 


aj{n) 


r(n + j + g) 

T{n-j + l)jl2R 


( 21 ) 
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and the remainder function Rn,i satisfies the bounds 

IP t W <■ ( \ ^ r(n + £+ ^) / \z\ 

' - \Tr\z\) |r(n-£+|)|£! 


i—n+^ 


cosh(9(z)) ( — 


Proof. We expand on the proof outlined in Watson [71 §7.3, p. 205]. The change of variables 
t = cos 0 turns the Poisson integral (13) into 


Jni^) — 


iz/2y 


r(n + l/2)r(l/2) J_, 


cos( 2 ;t)(l — t^Y dt . 


We split 


with 


Jn{z) = + Jn (z)), 


Jniz) = 


iz/ 2 y 


r(n + l/2)r(l/2) 




and J~ ifz) = Jfi {z). Now we change the contour in the integral for J+ towards a Il-shaped 
contour consisting of the line segment from —1 to — 1 + iN, followed by the line segment 
from — 1 + iN to 1 + iiV, and then the line segment from 1 + iN to 1. On that contour as 
well as on its convex hull, we have 3?(1 — t^) >0 with equality only at the end points of the 
contour. Hence we may choose the continuous branch of the square root function on the 
slit plane C\ (—oo, 0) which takes positive values on the positive real axis. For simplicity of 
notation, let us introduce the half-integer u := n — 1/2. The integral over the first segment 
of the contour becomes 


(2;/2)'^+V2 

r(u + i)r(i/ 2 ) 



e-^^e-^\2is + s^Yids 


{2'Kzyl‘^T{v -|- 1 ) 

gi(-2:+(v+l)7r/2) 

(27r2;)i/2r(u -|- 1) 



e-V(i + 
e-“u^(l - 




Yd 


where in the hrst identity we replaced zs by u and rotated the contour towards the line 
segments from 0 to zN, and used r(l/2) = In the second identity, we first pulled 

an integer power of i out of the integral and and then pulled half a power of i out of the 
integral without leaving the domain of definition of the chosen branch of the square root 
function. If 3?(z) > 0, then this last integral has a limit as N ^ oo. Similarly, the integral 
over the third line segment becomes 


(2:/2)'^+^/2 

r(u + i)r(i/ 2 ) 



Y^e-^Y-‘^is + s‘^Yids 


gi(2-(i/+l)7r/2 

(27rz)V2r(zy1) 


pzN 

i ' 


77/ 

e-Vil + ^ydu . 
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If ?ft(z) > 0, then the limit again exists. Moreover, the integral over the second line segment 
tends to 0 as —)■ oo. We therefore obtain 




r(p + i) 




iu 

Yz 


Y + e" 


(1 + 


iu 

Yz 


V 


du , 


where ojn = z — n7r/2 — 7r/4 = z — {v + and the contour has been changed to one 

along the real axis. 

The function (1 + auY is infinitely differentiable in u G [0, oo) for any a G C \ (—oo, 0]. 
Taylor’s expansion gives, for any ^ > 1, 


(1 + auY = 


t-i 

fc =0 


au}^ Y 

T{v + l-k) 


r(p + 1) {auy 


where r is the t'-th derivative of the function u i—)■ (1 + auY at some point uq G [0, u]. If 
i > u, then this derivative is proportional to a negative power of the function and thus 
attains its maximum at the point where the value of the function comes closest to the 
origin. This point equals u = —K(a)/|ap and the value of the I'-th derivative there equals 


r(p + i) 

T{v + l-l) V a J 


Hence we can write for the remainder term 



9{u) 


r(p + i) 

T{iy + l-£) 


a 



£l 


for some \0{u)\ < 1 that also depends on a. Note that 


3?(a) I |a —3ft(a)| |3'(a)| 

a |a| |q;| ’ 

and the latter equals 'ik{z)/\z\ if a = ^i/{2z) = ^iz/{2\z\Y- For > 0, we have that 
1 ± ^ is not real, hence we may insert the Taylor expansion into the integral for J^, and 
obtain for {z) the expression 


r-i 


{2ttz) 


1 


Ai=0 


T{u + 1 - k)k\ Jq 




'n f _I 

\2zJ^ Uz 


du + R{z) 


= {2ttz) 


-1/2 ^ r(i/ + fc +1) 


fc =0 


T{u + l-k)kl 


S) 


2z 


+ R{z) , 


where the remainder term R = Rn/ satisfies 

2 \i/2 r{u + £ + l) 


\Rn,eY)\ < 


£-u 


cosh(/5(z)) { — 


.7r|z|/ \TY + l-£)\£\ \ifi{z) 

We split the summation into even integers and odd integers to obtain 



ESTIMATES FOR CERTAIN INTEGRALS OF PRODUCTS OF SIX BESSEL FUNCTIONS. 


11 




k even 


- isf - 

k odd 

This completes the proof of the lemma. □ 

We specify our findings into explicit first order asymptotics valid for sufficiently large z 
near the positive real axis. We have with the notation from the above proof: 

Corollary 7. Let z € C be such that Q{z) < ?R-{z). Then 


/ 2 \ 1/2 1 / 2 \ f 

-[-zj - TI (^) ( 

/ 2 \ 1/2 o / 9 \ ^/^ / 

J cos(a;i) < — cosh(3=(2;)) ^ 


|3?(z)| 


|3?(z)| 


and if n > 1 and > n‘^, then 


Jni^)- ' — cos(a;n) < ' -p-r —r cosh(3'(2;)) 


1/2 


IK(z)l 


Proof. Following the previous argument, only the last inequality requires explanation. We 
choose i = n = n + 1/2, and apply the above expansion with the observation that 

Ffu + l + Zc) / l\2fc 

rUi-/ d-+2) 


for k < V. It follows that 


/ 2 \ 1^^ / 2 \ t _i 

cos{un) ^ cosh(3'(2;))(2|z|)“^ + i? 


cosh(/5(z))J^(l-2i ^) + R< 


cosh(3'(2:))^ 


/ z 


\Tr\z\J l^l KT^lzU kl \\^iz)\J 

where in the last line we estimated the remainder similarly to the terms in the expansion. 


One proves explicit asymptotics with one extra term in a similar way: 
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Corollary 8. Let z € C be such that 3f?(z) >1/4 and 9‘(-z) < Then 


Joi^) - 


2 \ 1/2 1 / 2 \ 

— cos(a;o) + — — sin(u;o) 

'KZ J az \TTZ J 


< 


< 


9 


128|zP VtIz 


1/2 


cosh(|3‘(z)|) 




5/2 


j±(,) _ ( A) ' eos(c.,) + 4 (A) ' sind^i) 


< 


< 


15 


128|z|2 Vvr|z 


1/2 


cosh(/j( 2 ;)) 


|K(z)| 


3/2 


and if n > 1 and 3?(z) > n^, then 


T±, . / 2 . 4n2-l / 2 y/2 

AL)-(-j cosK) + ^^^^-j smK) 


< 


1 / 2 

< - I I cosh(9=(z))— 


4 / II \ 1/2 
n / z ' ' 


4 \^7r|z|y ' ^ \l^(2)l 

Finally, we need zero order upper bounds for the asymptotic expansion: 

Corollary 9. For r > 0, we have that 
(a) 


9 / 2 y 


(b) 






Proof. This follows for r > 1 from Corollary while for r < 1 it follows from the trivial 
bound \Jn\ <1. □ 


Remark. Using more refined oscillatory techniques related to Sturm’s comparision princi¬ 
ple, the sharper bound Jo(t)| < (2/7r)^^^ is established in [3]. However, the bounds 
given by Corollary suffice for our purpose, and its proof is more in light with the elemen¬ 
tary nature of the present paper. 
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3. Useful estimates involving the Gamma function 
A version [3] of Stirling’s formula for the Gamma function states that, for x > 0, 

(22) r(x) = 

where the function ^ satisfies the double inequality 

(23) 




The starting point for all the convex estimates we need is the following well-known result: 
Lemma 10. For x > 0, the function x i—)■ log(r(x)) is convex. 

Proof. Let x,y > 0 and 0 < A < 1. Set P = j and q = It suffices to show that 


(24) 


r(- + ^) <r(x)rr(y)^, 

\p qJ 


for then the result follows by taking logarithms on both sides. To verify (24), consider the 
auxiliary functions 

31 — 1 _t_ 1 / —1 

f{t) := t p e p and g{t) := t p e <i, 

which satisfy 

fit)g{t) = f{tr = e-^e-\ and g{tY = ty-^e-\ 

The result is thus a consequence of Holder’s inequality. 

Corollary 11. Let x > Then: 

v(x + \^ 


□ 


y/x- 1/2 < 


r(x) 


< y/x. 


Proof. Use the identity x = r(x -I- l)/r(x) together with log convexity of T. □ 

Corollary 12. Let x > y > 0 and w >0. Then: 

r(x) ^ r(x -|- w) 
r(y) “ r(y-hu;)' 

Proof. For fixed tc > 0, we claim that the function 

r(x) 

r(x -i- w) 

is decreasing in x. This happens if the inequality 

r'(x) ^ r'(x -|- w) 
r(x) “ r(x -I- w) 

holds for every x > 0, which in turn is a consequence of the log convexity of T proved in 

Lemma [TOl d 
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The following lemma estimates binomial coefficients. 
Lemma 13. Let x > 1 and let 0 < d < x be an integer. Then: 


(25) 


r( 2 x) 


< 


el/24 


r(x — d)T{x + d) 2y^ 


,ll2^x^-(Tlx 


Proof. By Stirling’s formula (22), we have that 

r(2x) ^ u(2x)-2u(x) ^ f^^^^u{2x)-2^{x) 

r(x)2 2^ 

Also, since x > 1, we have that /i(2x) < and therefore 

^fi{2x)-2ti{x) ^ ^f^[2x) ^ 

It follows that 


\/2^2x 

r(x)2 “ 201 


By induction, observe that 

r(x)r(x +1) 


r(x — (i)r(x + d +1) + j 


^ d ^ TT ..-2ilx _ ^-{<P+d)lx 


j=i " i=i 

whenever d is an integer satisfying 0 < d < x. This is a consequence of the elementary 
inequality 

< e-. 

1 I — ^ ’ 

1 + 2 / 

which is valid in particular for y G [0,1]. Using (x + d)/d < we obtain that 


r(x)r(x) 


r(x - d)r(x + d) 

again for integers d G [0, x). Finally, 


< e^!^ -—- < e^!^ JJ e 

i=i i=i 


- 2 ilx _ g-+/i 


r(2x) 


r(2x) r(x)r(x) 


,1/24 


< 


r(x —d)r(x + d) r(x)2 r(x — d)r(x + d) 2 ^ 

as desired. 

Remark. The inequality 

r(x)r(x) 


^l/222x^-d?/x 


□ 


r(x — d)r(x + d) 


< e 


-tP/x 


is still valid for non-integer values of d G [0,x), and therefore Lemma 13 holds in this case 
as well. 
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We will also need estimates for the magnitude of the coefficients aj{m) defined in (21) 


when j is close to m. More precisely, in Section we will need good upper bounds for the 
quantities \am+A{fn)\. For small values of m, we compute them explicitly. For large values 
of m, we estimate these quantities via the following lemma. 

Lemma 14. For any natural number m > 1, we have that 

\am+i{m)\ < 

16 V TT 

Proof. The goal is to bound the absolute value of 

F(2m + 9/2) 


Rm+4(^) — 


F(-7/2)r(m + 5)2»"+4' 

Making repeated use of the identity r(x + l) = xr(x), together with the convexity estimate 
from Corollary 11 and Stirling’s formula ( 22 ), we have that 


\am+A{m)\ < 2 


5 F(2m-l/2) 

F(-7/2)r(m)2™+4 

1/2 r( 2 m) 


< 2(2m- 1) 


F(-7/2)r(m)2™ 
2(2m - 1 )-V 2 ( 2 m) 2 m-l /2 


m 


m—ll2 


r(-7/2)2"' 

The bounds (23) for the function /i imply that ^(2m) < and so 

follows that 


(26) 

The value r(— 7/2) 
completes the proof. 


('mil 2(2m 1) ! r,2m-l/2 m -m 

|am+4(m)| < ^ ' 

can be computed via the second formula in (14), and this 


□ 


4. Part I. Expanding Jq and Ji 

In the next four sections, we will be working under the standing assumption that n > 
no := 20. We start by asymptotically expanding the Bessel functions of order 0 and 1 and 
their relevant products. Due to need of accuracy, we must consider asymptotic expansions 
of length six and keep track of all the terms. The following notation will be convenient. 
Let us say that 

a ~ (ao) + (ai) + ( 02 ) + ( 03 ) + ( 04 ) + ( 05 ) with remainders fq, n, r 2 , ra, r 4 , rs, re 
if 
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a — 


k 


< Tfc, for every 0 < fc < 6 . 


2=1 

We also call rg the last remainder. Suppose additionally that 

b ~ ( 6 o) + (bi) + ( 62 ) + ( 63 ) + ( 64 ) + ( 65 ) with remainders sq, si, • 52 , S 3 , S 4 , S 5 , sg. 
Then we have the following product formula: 


with remainders 


5 k 

ab ~ '^C^aibk_i) 
k =0 i =0 


k 

rosk + E ri\hk-i\, for every 0 < /c < 6 . 
2=1 

Recall the coefficients ao(n) through 05 ( 71 ) for n = 0, 


9 

75 

3675 

59535 

2d012d5 

57972915 

128’ 

102d’ 

32768’ 

2621dd’ 

dl9d30d’ 

3355dd32 


and for n = 1 , 


15 

105 

d725 

72765 

2837835 

66891825 

128’ 

102d’ 

32768’ 

2621dd’ 

dl9d30d’ 

3355dd32 


To make the forthcoming notation less cumbersome, let us define, in view of Lemma 

^n{r) := (yj Jnir). 

To avoid writing many fractions, we further define t := (16r)“^. We also set 

c := cos(r — vr/d) and s := sin(r — vr/d). 

From Lemma and Corollary we have that 

ao(r’) ~ (c) + (2te) + (-ISt^c) + (-Booths) + + (2381d0t®s) 

with remainders discussed there as 

I, 2t, 18^^ 300^^ 7350t'‘, 2381d0^^ 960d980t®. 

8 

In a similar way, 

ai(r) ~ (s) + (etc) + (30t2s) + (-d20t^c) + (-9d50tS) + (291060t'^c) 


—, 6t, 30^^ d20^^ 9d50t'^, 291060^^ 113513d0t®. 
8 


with remainders 
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Applying the product formula, we obtain successively 

5o(t) ~ (c^) + (4tcs) + (—+ (—672t^cs)+ 

+ (15024t^c2 - 1200 tV) + (516480t®cs) 


with remainders 


81 17 169 2 1419 o 68571 . 1092495 . 43435485 . 

—, —t, —r, -r, -r, -r, -r, 

64’4’4’2’4 ’ 2 ’ 2 


(27) 50 ( 7 -) ~ 2ooo{r) :=(c^) + (6tc^s) + (-54f^c^ + 12t^cs^) + (-1116t^c^s + 8t^s^) 

+ (23022t^c^ - 3816t^cs2) + (836964t®c2s - 3600tV) 

with remainders 

729 217 2353 2 20003 , 956787 4 15017799 . 588969477 . 

—, — t, - r, - r, -r, - r, -r. 

512’ 32 ’ 32 ’ 16 ’ 32 ’ 16 ’ 16 


In particular, 

(28) |do(r) -dooo(r)| < - — -(16r) ^ 

On the other hand, we have 

5 i(r) ~ (s^) + ( 12 tcs) + (36t^c^ + 60t^s^) + (—480t^cs) + 

+ (-5040t^c2 - 18000tV) + (443520t®cs) 


with remainders 


121 57 429 2 2715 o 113535 4 1659735 . 62391105 . 

—, —t —r, -r, -r, --r, 

64’4’4’2’4 ’ 2 ’ 2 


('3i5o)(7’) ~ ■3iio(7’) ;=(cs^) + (12tc^s + 2te^) + (36t^c^ + 66t^cs^) + (—624t^c^s — ISOt^s^) 
+ (-5688t^c3 - 16290t^cs2) + (519480t^c2s + 184140tV) 

with remainders 

1089 577 5433 2 38331 , 1638411 4 23971455 . 897834285 . 

-, -1 , -1 ^, -- f. 

512 ’ 32 ’ 32 ’ 16 ’ 32 ’ 16 ’ 16 


In particular. 


^ ^ .M 897834285, 

|(d?do)(r) -diio(r)| < - — -(16r) ^ 


( 30 ) 
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Inequalities (28) and (30) are at the core of the following result, the proof of which does 


not require m to be even, nor m <n. 

Estimate A. For n > no = 20 and m > 0, we have 

roo 

(31) / Jn+mJn^m{^l - ^ooo)r~^dr < O.Tdng + m)“®, 


(32) 


I 

O 1 1 / 2 ^ 

Jn+mJndm{dlO0 illlo)l' dv < 1.12?T,q (tI -|- ?7i) 


Proof of Estimate A. Using the Cauchy-Schwarz inequality and Lemmataandto com¬ 
pute the integrals, we have that 


1 \ 1/2 /128 r(n -h m - ^) \ 1/2 

i) \693 rln-I- 777,-L hli / 


dr \ 1/2 


_ / 1 \ 1/2 / 128 r(n -h m - ^)' 
\2n/ V693 r(7r-|-771-|-^)^ 


where 


r(£-^) ,2 r(f) ,2 

ao := sup —-= , -V 20^^ < 1.21. 

£>nor(^ + f) r(f) - 

Thus, for n > no, the left-hand side of inequality © is bounded by 


where 


588%9477/16 y.. 64y/2 
166 0 V693/ 


= -;;-an - ^ U.l^. 

16® ° V693/ 

In a similar way, for n>no, the left-hand side of inequality (32) is bounded by 


where 


^ 897834285/16 1 / 2 /64 n 1/2 

^ =-“» (™1 


< 1 . 12 . 


□ 
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5. Part 11. Expanding Jm 

Start by noting that cosum = (—and sincj^ = (—since m is an even 
integer. We have that 

(33) 


f+i 


-+i 


- Pmir) = (-l)'^(16t)2%,(m) -s (-l)'=(16t)2"+'a2fc+i(m)), 

k=0 k=0 

where the error term pm is implicitly defined by this identity. From Lemma we know 
that 

\am+A{m)\ 


(34) 


\Pm{r)\ < 


r-m+4 


for every r > 0. The main integrals, which are the subject of the next section, arise from 
replacing HmUooo by 


(35) 

and HmUiio by 


■vlmOOO •— (dm Pm)il0005 


(36) 


dmllO •— (dm Pm)dll0- 


Estimating the error term in this replacement is the main goal of the present section. We 
formulate it as 

Estimate B. Let n > no = 20 and m be even. 

If m = 0, then 


(37) 

(38) 

If m = 2, then 

(39) 

(40) 

If m = A, then 


poo 

/ <7n(5o5ooo — 3 oooo)t ^dr 

Jo 

poo 

/ «^nDo-3llO — 3oilo)T“^dr 

Jo 


< 0.022no ^n"^. 


< 0.023nQ ^n 


poo 

/ Jn+ 2 Jn {^22000 — ^ 2000 )^ 

Jo 

poo 

/ Jn+2Jni^2^nO — ^2llo)'l^~^dr 

Jo 

poo 

/ <7n+4<7n(345oOO — 5400 o)f dr 

Jo 


< 0.162no ^n“®. 


<0.166nQ^n 


< 2.823nQ ^n 
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Jn+iJniZiZiio — Ziiio)r ^dr 


< 2.885no ^n"®, 


If 6 < m < n, then 
(41) 


(42) 


poo 

I Jn+mJn{-iimd000 ■JmOOo)^ dr 
Jo 

POO 

I Jn+mJnidmdllO ■Jmlio)!’ dv 


< 0.015no ^n"^ 


< 0.015no ^n“^. 


Proof of Estimate B. We consider the case m = 0 first. The left-hand side of (37) is 
bounded by 

poo 

(43) |a4(0)| / J2(r)(l + 6t + 66t2 + 1124t®-h26838t^ + 840564t®)r“'^dr, 

Jo 

whereas the left-hand side of (pSl) is bounded by 


(44) |a4(0)| / J^{r){l + Ut + 102t‘^+ 80At^+ 21978t‘^+ 703620t^)r-^dr. 

Jo 

For i G {5,6,..., 10}, Lemma implies that 


2 -^r(£)r(n-^) 
/ J„{r)r dr = — y —— 
lo r(^)2r(n + ^) 


This can be estimated as follows: 


poo 

j J^{r)r~^dr < 


where 


( 0 ) 2-^m r(20-^) 


cl := 


r( 


t +1 


2 ;^r( 20 + 2 


£+ 1 - 


20 ^ 


In particular, for £ G {5, 6,..., 10}, one can easily check that 


0.14 < c® 


< 0.19. 

It follows that the upper bounds for each of the last five summands on the right-hand side 


of inequalities (43) and (44) can be estimated by a small fraction of the upper bound for 


the first summand. Quantifying this, one obtains (37) and (38), respectively. 

We consider the case m = 2 next. Again for i G {5, 6 ,..., 10}, the integral to consider 
is the following: 
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\Jn{r)Jn+ 2 {r)\r ^ 


In view of the absolute value in the integrand, this cannot be computed directly with 
Lemma 1^ Instead, we use the Cauchy-Schwarz inequality to estimate 


f 


\JnJn+2\r ^ ^dr < 


“ 2 dr y/2 

'^n+2^2e+3 ) 


0 ^ ^ ^^0 

1 \ 1/2 ^2-(2^+3)r(2£ + 3) r(n - ^ + 1)1/2 


r(^ + 2)2 r(n + ^ + 4). 

where the last identity is a consequence of Lemmata and Reasoning as before. 


we 


derive the estimate 


\JnJn+2\r ^ ^dr < cf^n ^ ^ 


^(2) / 2-(2^+3)r(2£ + 3)r(20-£ + i) ^^2£+3 


where 


2 r(£ + 2)2 r(20 + £ + 4) 

In particular, for every £ G {5, 6 ,..., 10}, one checks that 


1/2 


0.11 < cf^ < 0.15. 

As before, it follows that the last five summands on the right-hand sides of estimates 


roo 


Jn+ 2 Jr 


i(-323ooo — 52000 ) 1 " 


^dr 


< 


and 


5n-|-25n(525ll0 — 5211 o)f ^dr 


< 


I»6(2)|(< 




14 




( 2).„-8 


102 


n 




(2).„-9 


804 


n 




(2).„-10 


n 


21978 


(2).„-ll 


104 


-9 


n 


+ - 


703620 

16^5 


( 2).„-12 


MO 


n 


can be bounded by a small fraction of the first summand. Quantifying this yields (39) and 

(4§. 

The cases m = 4,6,8,10 can be treated in a completely analogous way to what was done 
for m = 2. We omit the details, but remark that for m = 6,8 ,10 this method produces 
estimates which are stronger than (41) and (42). However, the latter will be enough for 
our purposes. 
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Finally, we deal with the case of even m > 12. Again for i G {5, 6 ,..., 10}, the integral 
to consider is bounded by 


roo 

i r(m + £)2 


1 ^^l/2^2-(2”"+2^-^)^(2m + 2£-1) r(n-£ + l) \V2 

r(n + 2 m + f')/ 


Identifying a binomial coefficient, we notice the trivial bound 

2-(2m+2£-i)p(2m + 2^-1) ^ 1 
r(m + £)2 ' - 2' 


It follows that 


-^-^dr < r(n-.g + l) \V2^ 

vr(n + 2 m + £) / 


\JnJn+m\'^ dv ^ Tl 


To handle the coefficients |am+ 4 (^)|, we recall Lemma 14 For even m > 12, define 

' ■ 2V7r^ ^ ln2£j-^2(20-^ + l + A:) 

a decreasing function of m for fixed L We finally arrive at 


1/2 


f 


dn+mJn(iimxf000 ilmOOo)^ 


< 


105 / (m) —5 6 (m) —6 66 (m) —7 1124 (m.) _g 26828 (m) —9 840564 _ 


16 "^ + 161 
and 


Co m +- 


164 


-Cq 'n + 


165 


cio'n 


10 


f 




105 




14 




102 


16 "^+ 16 ^“' "^+ 16 i >‘=« "^+^^''<1 




< 


804 




21978 




Cq 'n +- 


703620 


165 


(m) -10 

c\o n 


If m > 12, then both of these expressions are bounded by 0.015nQ as claimed. □ 


6 . Part III. The core integrals 
The main integrals that are left to analyze. 


fJ,Q = fiQ{m,n) := / Jn(l’)Jn+m(T) 3 m 000 ( 7 ’)r dr] 
Jo 


fii = Hi{m,n) := / Jn{r)Jn+mir)^miio{r)r dr, 
Jo 
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decompose into the core integrals Q by means of expanding 3 m 000 and 3miio using the 
following elementary trigonometric facts: 

8 c^ = — cos(4r) + 4sin(2r) + 3, 

8 c^s = — sin(4r) — 2cos(2r), 

8 c^s^ = cos(4r) + 1, 

8 cs^ = sin(4r) — 2cos(2r), 

8 s^ = — cos(4r) — 4 sin(2r) + 3. 

These identities can be readily checked recalling the definitions c = cos(r — vr/4) and 
s = sin(r — vr/4), and noting again that cosujm = (—and sinwm = (— 

A simple parity check verifies that the resulting core integrals with cos(2r) and sin(2r) 
satisfy the parity assumption of Lemma relative to the powers of r, and so these terms 
yield zero contribution. It therefore suffices to consider the constant terms and the terms 
involving cos(4r) and sin(4r). The strategy will be to split the main integrals 


(cos) , (sin) 

Mo = pto + Mo 

(cos) , (sin) 

Mi=Mi +Mi 


according to cosine and sine contributions. More precisely, recall definitions (35) and (36) 
for 3 m000 aud 3miio, respectively. The first factor in each of them, namely — Pm, is 
given by identity (33). The right-hand side of this identity consists of two sums which 
come affected by a coefficient of c and s. These are at the source of what we denote by 
cosine and sine contributions, respectively. Working out the algebra, one is led to define 


(45) p 


(cos) 


^ (-l)'‘a2fc(m) / Jn{r)Jn+m{r) ■ {ao + a 2 t^ + a4t^)+ 


k=0 


+ (/3o -h cos(4r) (yit -h + 75 ^^) sin(4r) 




dr 


for the sequence of coefficients given by 


(46) 


(aoj 0:2, 0:4) 


(3,-150,65250), if =t= = 0, 
(1,174,-33354), if * = 1; 


(47) 


(/3o,7u/32,73,^4,75) = 


and 


(-1, -6, 66,1124, -26838, -840564), 
(1, -10, 30,444, -10602, -335340), 


if * = 0 
if * = 1 
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, . , /I 'jm/2 

(48) := -^(-l)^a2fc+i(m) / Jn(r) • (ait + ast^ + Q;5i®) + 

JQ L 


^+1 


8 


fc =0 


+ (/?it + /Sst^ + /Sst®) cos(4r) + (70 + 72 t^ + 74 t‘^) sin(4r) 


-2k-2 


dr 


for the sequence of coefficients given by 


(49) 


(ai, as, as) 


(6,-1092,826164), if * = 0, 
(18,-1164,1071900), if* = l; 


(50) 


(70,/3i,72,/?3,74,/35) 


(-1, 6 , 66 , -1124, -26838,840564), if * = 0, 
(1,10,30, -444, -10602, 335340), if * = 1. 


For * G {0,1}, the goal is to obtain a set of estimates of the form 


(51) - Mi“®^(m,n)| < E[^°^\m,n) + 

(52) (m, n) - n)\ < n) + n), 

where M and Ei denote the main and error terms coming from the analysis of the constant 
terms, and E 2 denotes the error term coming from the analysis of the terms of frequency 
4r. We shall denote Ei and E 2 by error terms of the first and second kind, respectively. 

6.1. Constant terms. Everything originating from the constant terms can be explicitly 
computed, one just needs to be careful about bookkeeping. As indicated before, we organize 
the terms into cosine and sine contributions. 


6.1.1. Cosine contributions. We split the analysis in four cases: m = 0, m = 2, m = 4 and 
m > 6 . In each of these four cases we will identify, as announced, a main term M and an 
error term Ei. 

Let us start by handling the case m = 0. In this case, the contribution coming from the 
non-oscillatory term ao + a 2 t^ + a 4 t^ in (45) can by computed exactly with the help of 
Lemmata]^ and 1^ the result being a main term 


Mi“"^( 0 ,n) 


1 1 

— an- “t“ 

8 2 n 


1 

8 


and an error term of the first kind 



1 

l)n(n + 1) ’ 


^(cos) 


(0,n) 


1 

8 


«4 a 2 \ 3r(n - 2) 1 

164 “2lUj^g2j i6r(n + 3) 8 


02 ( 0 ) 


a 4 \ 5r(n — 3) 
W) 32r(n + 4)' 
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Recalling (46), the main term can be computed as follows: 
(53) 


(0, n) = <j 2048(n-l^)n(n+l) > 


51 


, -*• I _ 

t 16n ^ 2048(n-l)n(n+l) ’ 

To estimate the error term, we first compute it as 


if * = 0, 
if * = 1. 


{ 101925r(n-2) 1468125r(n-3) -n _ p, 

4194304r(n+3) 1073741824r(n+4) ’ II * — U, 

54729r(n-2) 750465r(n-3) -r _ . 

4194304r(n+3) 1073741824r(n+4) ’ U * — i- 

Using the triangle inequality together with the easily verified bounds 

1 r(n-2) 1.02 , 1 r(n-3) 1.04 

n° 1 (n + 3j n‘ 1 (n + 4j n' 

valid for n > 20, one arrives at 


(54) 


1,* u - 1 o.015no^re 


-i„-4 

-1^-4 if* = l. 


We move on to the case m = 2. Orthogonality kicks in the form of Lemma to ensure 
that we only have one main term, which the same lemma computes as 


Mi“^)(2,n) := 1 ( - ao(2)^ + a2(2)ao) 


1 


/: 


195 


In other words, 

(55) _ _ 

4096n(n+l)(?iH-2)' 

The error term of the first kind is now given by 
1 


Mi“®^(2,n) = I 4096n(n+l)(n+2) 

I ‘ 


8n(re + l)(n + 2)' 

, if * = 0, 
if * = 1. 


r(n - 1) 


„) :=- ( - + ,..(2)^ - a.(2)„„) ^ 

^ “4 „ “2 \ 15r(n-2) 1 / ^ 04 \ 7r(n-3) 

+ 8(^®2(2)i 64 i28r(n + 5) 8 V 164/ 64r(n + 6)' 

Proceeding as in the case m = 0, we see that this term obeys the following estimate: 


(56) 


l^iT^(2,n)| < 


0.039no if * = 0, 

0.012no if * = 1. 


In the case m = 4, we expand to one higher order. The reason for this will become 
apparent at the end of Section We thus have exactly one main term 
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and an error term of the first kind 


1 

32n(n + i)(n + 2)(n + 3)(n + 4) ’ 


j^(cos) 


(4,n) 



+ 04(4:) 
- 06(4) 


02 \ 7r(n - 2) 1 / 

l^Jl28r(n + 7) 8 V ’ 


04 \ 15r(n — 3) 
256r(n + 8)' 


As before, we compute the main term 


r_ 322425 _ if* = n 

(57) ^ I 1048576n(n+l)(n+2)(n+3){n+4)’ 


_7(m_ 

L 1048576n(n+l)(n+2)(n+3)(n+4) ’ 

and verify the following bounds for the error term: 


if * = 1 , 


(58) 


|i?fr^(4,n)| < 


- 1^-6 


0.42nQ n 


if * = 0 , 


O-IIrq if * = 1 . 


If m > 6 , orthogonality ensures that there is no main term. The error term of the hrst 
kind is given by 


r(n) 


2m+ir(n + m + 1) 

(m + 2)r(n — 1) 


T-,(cos) / \ I / / \ 0^4 / \ 0^2 . / 

-^ 1 ,* [m,n) ■=^\am-A{m)r^ - 0 ^- 2 ( 711 )^ + amijn 

If , ^ 04 / N «2 / N \ (m + 2 )i(n-i] 

1 / 04 , , 02 \ (m + 3)(m + 4)r(n - 2) 

+ g(o™(7o)^ - am+2{m)^) 2-+6r(n + m + 3) 

If 04 (m + 4)(m + 5)(m + 6 )r(n — 3) 

+ 3 .2™+8r(n + m + 4) ’ 

and this can be crudely bounded in the following way. If m > 6 , then 


(59) 


^fr^(m,n)| < |i?17'(6,n)| 


(cos) , 


for the given range of admissible m and n. It is easy to see that inequality (59) holds if 


m is large enough, essentially because each of the summands that constitute the left-hand 
side of that inequality is of order at most For the remaining cases, one checks it 

directly. The upshot is a bound of the form 
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(60) 


l^(cos) 


(m, n)| < 


6.34ng 

O.OOng 


if * = 0, 
if * = 1, 


valid for every even m > 6. 


6.1.2. Sine contributions. We proceed similarly, again splitting the analysis into four cases. 
If m = 0, then the contribution coming from ait + a^t^ + a^t^ amounts to a main term 


M, 


(sin) 


1 


(0,n) := --ai(O) 


ai 


1 


8 16 4(n — l)n(n + 1) ’ 


and an error term of the first kind 


^ > SV ^163 ^16/16r(n + 3 


+ 


8 


(ai(0) 


as 

16^ 


16/ 16r(n + 3) 

"3 5r(n - 3) 1 

'“^^^^163]32r(n + 4) 8 


Recalling (49), we compute the main term as 


as \ 35r(n — 4) 
lf)3/ 256r(n + 5)' 


(61) 


Mi““y0,n) = { 2048(n-l)n(n+l) 


, if * = 0, 
if =t= = 1. 


f 2048(n-l)n(n+l) ’ 

Arguing as in the last subsection, the error term can be seen to obey the following bounds: 


(62) \E 


(sin ), 


0.0016no^n-'^, 


if * = 0, 


1 ,* ( 0 ,n)| < 'v „ „„„„ -1 _4 .p , 

lO.OOdOng n , if * = 1. 


If m = 2, there is a main term 


(63) Mf"y2,n):=^(ai(2)^V 


1 


/ 


45 


_ I 4096n(n+l)(n+2) 


, if * = 0, 


167 8n(n + l)(n + 2) \ 4096n(nTl)(n+2) > = 

and an error term of the first kind 


-E 


(sin) 

1 ,+ 


fo 1 701 “3 , ai\ r(ri-jO 

(2,n). + ^^^^6/8r(n + 4) 


701 "5 „ «3 ai\15r(ri-^ 

'^8V ^^^^165 ^*'^^163 ®^^^16/128r(n + 5) 

1 

+ A 


I li _ 

165 ^163;64r(n + 6) s' 


V ®^^165/ 1024r(n + 7) 
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which satisfies 

(64) |i?i7^(2,n)| < 

If m = 4, we again expand to one higher order. There is a main term 


'o.0062no^n-^, if * = 0, 
0.0031no^n-^, if * = 1. 


“^^^^16/32n(n + l)(n + 2)(n + 3)(n + 4) 

and an error term of the first kind 


1 


as 


.oiy 


1 


^(sin) 


(4,n) 


ai \ 3r(n — 1) 


+ ^ (- + “= Wll - 

y (<..{4)(| - y ( - a,(4) 


as \ 495r(n — 4) 

165j8192r(n + 9)' 


The main term satisfies 


(65) 


_ 76167 _ 

1048576n(n+l){n+2){n+3){n+4) ’ 

_ 211869 _ 

1048576n(n+l){n+2){n+3){n+4) ’ 


and the error term can be bounded as follows; 


if * = 0, 
if * = 1, 


( 66 ) 


d:;"’(4.»)i < 


0.086nQ 

0.063no^n-^ 


if * = 0, 
if * = 1. 


Finally, if m > 6, there is no main term, and the error term of the first kind is given by 


EiT’ 


r(n) 


V ’ ^ 8V ^65 ^103 ™ ^16/2”4+ir(n + m + l) 

1/ . X as / ^ “3 , / xai\ (m + 2)r(n-l) 

1/ / ^ “5 , / N “3 / ai\ (m + 3)(m + 4)r(n - 2) 

+ 8i iM^gS+«-+iM^g3 «m+3(m)^J 2™+6r(n + m + 3 ) 


/ N “5 , as \ (m + 4)(m + 5)(m + 6)r(n - 3) 

+ o ( ®»Tt+l(7n.) p. am+3(7n) .. „Q ) „ om.-l-S-n('„ I ™ 1 


1 


165 


163 7 3 • 2™+8r(n + m + 4) 


1^ as (m + 5)(m + 6)(m + 7)(m + 8)r(n — 4) 

+ -( -a,„+3(m)^j 3 • 2”4+i2r(n + m + 5) ' 


Again the monotonicity formula 
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< l<“’(6,n)| 

holds for every even m > 6 , and this implies a bound of the form 

,- 3^-4 


(67) |Ej7)(m,n)|< 

which is valid in that range of m. 


1.49nQ'^n 

4.08ng 


if * = 0 , 
if * = 1 , 


6.2. Frequency 4r terms. To handle the terms of frequency 4r, we make repeated use 
of the following result: 

Proposition 15. Let n, m G N he such that n > uq = 20 and m even with 0 < m < n. Let 
a G {1,3,5} and jd G {2,4,6}. Then each of the following quantities is less than n“^0.35”'.‘ 


(*) 

{ii) 

{Hi) 

{iv) 


—+1 

2 poo 

V / Jn(r)Jn+m(T)r“^*^a2fc(m)cos(4r)r“"(ir 
Jo 


fc =0 

—+1 
2 

E 

—+1 

E 


Jn{r)Jn+m{r)r ^^a 2 fc(m) sin(4r)r ^ dr 


Jn{r)Jn+ra{r)r ^a 2 fc+i(m) cos(4r)r ^ dr 


fc=o 
—+1 

2 /*00 


^_ 1*00 

^ / Jn(r)Jn+m(r)r“^^“^a2fc+i(m.)sin(4r)r““dr 

fc=o 


Proof. All estimates can be proved in a very similar way. We focus on the tightest case, 
that of (z) with a = 1 , and briefly comment on the other cases at the end of the proof. 
Using the definition of the coefficients aj{n) and Lemma together with the convexity 
estimate from Corollary 12, we obtain 


^ |a 2 fc(m)| / Jn{r)Jn+mir)r~‘^’" cos{Ar)r~^ dr 


k=0 


< 


< 


< 


E 

k=0 

—+ 1 
2 

E 

k=0 


T{m +2k+ 1/2) 


•,2k 


r(2n + m — 2k) 


T{m -2k + l/2)22fcr(2A: + 1) 42 ^+”^ r(n + l)r(n + m + 1) 


T{m +2k+ 1/2) 


r(2n + m) 


i—2n—m 


T{m — 2k + l/2)r(4A: + 1) r(n + l)r(n + m + 1) 


( 68 ) 
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The second fraction in this expression resembles a binomial coefficient and does not depend 
on k. It can be estimated in the following way: using Lemma 13 with x = n + m/2 and 
d = m/2, we see that 


r(2n + m) 


1 


r(2n + m) 


r(n + l)r(n + m + 1) n(n + m) r(n)r(n + m) 

1 eV24 


< 


(69) 


< 


n(n + m) 2^/^^ 

1 

2^/^ n(n + m)^/2 


(n + m/2)i/222n+mg-^^^ 


;)2n+m 


To estimate the sum of the first fractions in (68), we proceed as follows. For k = m/2 + 1, 
we simply have that 


r(m + 2fc + l/2) 


r(2m + 5/2) 


< 


r(m -2k + l/2)r(4/c + 1) r(-3/2)r(2m + 5) 

1 r(2m + 3) 


< 


r(-3/2) r(2m + 5) 4Vvr (2m + 4)(2m + 3)' 


On the other hand, as long a^O <k< m/2, we can use Corollary 12 followed by Lemma 
13 with X = m/2 + fc + 1 and d = |3A: — m/2| to conclude that: 


E 


r(m + 2A: + 1/2) 


< 


^ r(m + 2fc + l) 


E 


^_g r(m — 2 k + l/2)r(4fc + 1 ) r(m — 2k + l)r(4A; + 1) 

m 

r(m + 2fc + 2) 


1 


k=0 


r(m — 2k + l)r(4A; + 1) m + 2fc + 1 


< ^ (m/2 + A: + l)V^ ^fc^-(|L^ 


vr m + 2A: + 1 

fc =0 


For 0<A:<m/2, itis easy to check that the quantity ^ decreasing in k. 

Moreover, for A; = 0, we have that 

m + 1 “ ^ ^ 

Using this, we are left to estimate the Gaussian sum 


^This does not work for k = m/2 + 1 because the assumptions of Corollary 
conclusion fails. 


12 


are not met and the 



























ESTIMATES FOR CERTAIN INTEGRALS OF PRODUCTS OF SIX BESSEL FUNCTIONS. 


T 


m 


We start with the trivial estimate 


m 

±4^e- 

k=0 


(3fc-m/2)^ 
Tn/2 + /c + l 


T 


m ^ 


2 

E 

fc =0 


, (■3k-m/2y 

4 ^g m+1 


Changing variables of summation ‘il = ?>k — we see that 

£eL 

where L is the new summation set, given by 

(mm m'y 1 ^ 

L := { -,-h 1,..., — > C -Z. 

I 6 ’ 6 ’ ’ 3 / 3 

We estimate this sum by the product of the largest term and the number of terms 
m/2 + 1. To detect the largest term, define the function 

_(3W 

ipm{x) ■■= 4^e "*+1 . 

The unique solution xq G y] to the stationary condition ip'^{xo) = 0 is given by 

log 2 

xo = —^(m + 1), 

for which we have 

+^(xo)=4l-+i, 

where 

^ log 2 /log 2x2 

yl := 4 9 e 1 3 ) < 1.06. 

Since +m(^) < +m(xo) for every i £ L, we thus have that 

T™< + 

It follows that 


-+i 


( 70 ) 


E 

k=0 


r(m + 2/c + 1/2) 
r(m — 2k + l/2)r(4A: + 1) 


< 


< 


2ei/24 


< 


/TT 

103 2e^/24 


100 


(m + + l) ^ 


4^/^^ (2m + 4) (2m + 3) 
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where the last inequality holds since the second summand on the second line amounts to 
at most of the first summand. Finally, estimates (69) and (70) together imply that 


"o" H” i 

^ poo 

'^\a2k{m)\ / Jn(r)Jn+m(r)r“^^cos(4r)r“^ dr 

k=o 


-+i 


< 


< 


E 

fc =0 

103 26^/24 


r(m + 2/c + l/2) 


T{m — 2k + l/2)r(4/c + 1) / Vr(n + l)r(n + m + 1) 


r(2n + m) 


(m + 1 )-V 22 ™( 


=1/24 


—2n—m 


.100 0r 
< n-^2^/^A^2-‘^^. 

Since m < n, we hnally get the desired estimate: 


2 y^ n(n + nipl‘^ 


2n-\-m\ A—2n—7n 


< n-i2"/3^’"2-2" = < n-^0.35". 

This completes the estimate of sum (/) with a = 1. For the other cases, letting k G 
{0,1,..., m/2 + 1} and 1 < j < 7, one just checks that the bounds given by Lemma 
namely 


poo 

/ Jn(r)J„+m(r)r“2^“-^e^*'’dr 

Jo 


< 


22 k+j-i ^2n + m-2k-jy. 


42 n+m n!(n + m)! 
are decreasing in j as long as the conditions of the statement are met. 


□ 


Remark. We will need the following observation for the purpose of our applications. For 
n > 20, we have that 0.35"'/^ < n~^, and so the bound given by Proposition 15 can be 
further estimated as follows: 


n“i0.35" = re“^0.35"/^0.35’"/2 < 0.35"«/^n“^ < 0.6”°n“^, 

provided n > no = 20. Alternatively, still for n > 20, we have that 0.35'^"' < n~^ if 
r > 0.72. Using this bound instead, we see that 


n-iO.35’" = n“^0.35™0.35(^“^)’" < 0.35(^“^)’"°n“® < 0.75'^On“®. 

All in all, we have the following upper bound for the quantities considered in Proposition 

dsi 

^ min{ 0 . 6 ”°, 0.75"'°n“^}^ n~^. 

This distinction will play a role to ensure good bounds for the m = 4 terms which were 
expanded to one higher order in the last subsection. 


We are finally ready to estimate the contribution coming from the oscillatory terms 
(/3o + /d 2 ^^ + cos(4r) and (yit + 73 /^ + 75 /^) sin(4r) in expression (45), and similarly 
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in (48). Appealing to Proposition 15 and the remark following it, we see that we can take 
the following for errors of the second kind: 


17 (cos) 
2 ,* 



16 162 ^ 163 ^ 164 ^ 165 A 




IM I M I IM I M 

16 162 163 164 





where {6,t) = (0.75,6) if m = 4 and {9,t) = 
from (47) and (50), we obtain the estimates 


(0.6,4) if m 7 ^ 4. Plugging the values of /3 ,7 


(71) 

and, if m 7 ^ 4, 

(72) 


47^(4, n)|,|i?^7(4,n)|< 


(sin), 


0.39 • 0.75’"“n 
0.30 • 0.75’"0n 


-6 

5 

-6 


(m,n)I < 


(sin) / 


0.39 • 0.6"'°n 
0.30 • O.O'^On 


-4 


-4 


if * = 0, 
if * = 1, 


if * = 0, 
if =t= = 1. 


7. Putting it all together 

In the last section we analyzed the core integrals, which were decomposed into main 
terms and error terms. We derived several estimates which are recalled below in each case. 
These are used together with Estimates A and B to yield appropriate bounds, which are 
then evaluated at no = 20: 


(i) If m = 0, then we use the knowledge about the main terms coming from (53) and 
(61), the estimates for the error terms of the first kind contained in (54) and (62), 


and the bounds for the error terms of the second kind from (72), to conclude that 

3 


roc o 1 


1 


-4 


/ J„(r)3n(r)r ^dr -h , . , , 

/o 16 n 128 (n-l)n(n +1) 

< (^(0.022 + 0.026 + 0.0016)no ^ + 0.74n“®/^ + 0.78 • 0.6’^°) n"^ < 0.0030n 
and that 

[ Jnir)^iir)do(.r)r-^dr - ~ . 

/o 16 n 128 (n — l)n(n + 1) 

< (^(0.023 + 0.015 + 0.0030)no ^ + 1.12no®/^ + 0.60 • 0.6’^°) n"^ < 0.0028n"^. 


(ii) If m = 2, then we use the knowledge about the main terms coming from (55) and 
(63), the estimates for the error terms of the first kind contained in (56) and (64), 
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f 


and the bounds for the error terms of the second kind from (|72|), to conclude that 

15 1 


Jn+2[r)J„{r)Mr)Ur)r- dr - 

< (^(0.039 + 0.0062)no ^ + 0.74n“^/^ + 0.162no^ + 0.78 • 0.6”°)< 0.0028n"'‘, 
and that 

9 1 

/ Jn+2{r)Jn{r)Z2{r)Z\{r)Zo{r)r~^dr - ■ 

Jq 256 n(n + l)(n + 2) 

< (^(0.012 + 0.0031)no ^ + 1.12no + 0 . 166710 ^ + 0.60 • O.O^o^n""^ < 0.0015n"^. 

(iii) If m = 4, then we use the knowledge about the main terms coming from (57) and 
(65), the estimates for the error terms of the first kind contained in (58) and ( 66 ), 
and the bounds for the error terms of the second kind from (|71[), to conclude that 


f 


1557 


Jn+4(r)Jn(r)34(r)5o(0^ ^ iqW 

4096 n(n + l)(n + 2j(n + 3)(n + 4j 

< (o.7in~^^‘^ + (0.42 + 0.086)no ^ + 2.823710^ + 0.78 • 0.75”°)n“® < 0.197n“®, 
and that 

[ Jn+4{'f’)Jn{'r’)24{r)2i{r)^Q{r)r~^dr ---—___—.— -- 

Jo 4096n(n + l)(n + 2)(n + 3)(n + 4) 

< {l.l2no^^‘^ + (0.11 + 0.063)no ^ + 2.885no^ + 0.60 • 0.75”°)n"^ < 0.26477“®. 

(iv) If 77T, > 6 is even, then there are no main terms, and we use the estimates for the 
error terms of the hrst kind contained in (60) and (67), and the bounds for the 
error terms of the second kind from (72), to conclude that 

poo 

/ Jn+mir)Jnir)^m{r)^l{r)r~^dr 

Jo 

< (o.015no ^ + 0.74n“^/^ + (6.34 + 1.49)770^ + 0.78 • 0.6”°) 77 “^ < 0.002277“^, 
and that 


Jn-\-m {r)Jn{r)^m{r)Zl{r)Zo{r)r ^dr\ 

< (^0.015no ^ + 1.12770 + (0-09 + 4.08)77o^ + 0.60 • 0.6”°) 77 “^ < 0.002077“^. 

To get the constants promised by Theorem one just multiplies the far right-hand sides 
of each inequality by the normalizing factor^ < This concludes the proof of Theorem 
[T]for77>20. 
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8. Numerical estimates for n < 20 

In this section, we numerically evaluate the integrals Iq and Ii defined in ([^ and Q, 
respectively, for 2 < n < 19 and even 0 < m < n. We split the integrals into 

rR roo 

Ij — Ij^\ow “ 1 “ -^jjhigh — / • • • dv I ... dv . 

Jo Jr 

We use a quadrature rule for the first integral and estimate the second integral by analytic 
methods. We aim at an absolute error of at most 0.9 x 10“® for Iq and Ii. 

The high integral would be entirely negligible at our desired accuracy for the threshold 
(say) R = 10^*^, but this would put unnecessarily much computing time on the low integrals. 
We choose R = 63000 and estimate the high integrals by a more careful analysis of the 
asymptotic expansion. To bring down the computing time for the low integrals, we use a 
high degree Newton-Coates quadrature rule. 

We first discuss the high integrals and begin with /o,high- Since R is large compared to 
(n + m)^, we take advantage of the asymptotic information in Corollary]^ Splitting each 
Bessel function into main term plus error, and applying the distributive law, yields one 
main integral of the form 

3 


main,high 


fR 


TTr 


- cos ( Wn+m ) COS {uin ) COS {Um ) COS'^ (wq ) vdr 


plus 2® — 1 error terms. 

If n is even, since m is even as well, an even number of the integers n,m,n + m is 
congruent two modulo four, and we obtain, with the periodicity cos(wn) = — cos(a;n+ 2 ), 


/n 


0,main,high — .^0,main,high,even 


nCK. 

Jr 


2 y f. 

— cos^((jJo)rdr , 
7rr J 


a term which is in fact independent of the particular even n and m. If n is odd, then we 
obtain similarly 


ho,main,high — ho^main,high,odd 

Likewise, if n is even, we have 


hi,main,high hi main,high,even 


and, if n is odd, 

hi, main,high — hi main,high,odd ■ — 


IR 


IR 


IR 


-) 

TTr ) 


— ) cos^(a;i) cos^(a;o)rdr . 


-) 

Trr ) 


) cos^(a;i) cos^(t<;o)Tdr , 


2 y 

— ) cos'^(a;i) cos^(a;o)rdr . 

7rr / 


These integrals have closed-form expressions in terms of trigonometric and trigonomet¬ 
ric integral functions. Mathematica calculates these expression and evaluates them with 
prescribed accuracy, resulting in 

I ho,main,high,even 


- 1.2798 X 10“®| < 10“^°, 
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1-^0,main,high,odd “ 0.2560 X 10 < 10 

|A,main,high ” 0.2560 X lO-^] < IQ-^O, 

where the distinction between even and odd n is not visible at the prescribed accuracy in 
the case of /i,main,high- A Sample Mathematica code used to evaluate /o,main,high,even is the 
following: 

N[Integrate[(2/Pi)^ * Cos[r — Pi/4]® * r^“^\ {r, 63000, Infinity}], 20] . 

We now estimate the 2® — 1 error terms of /i,high — A,main,high- Of these error terms, 
six of them consist of an integral of a product of five main terms of Corollary and one 
error term of Corollary To estimate these six terms, we use the finer information from 
Corollary for the error term of Corollary 

The second main term of Corollary leads to integrals of the type 

^ sin(a;„+m) cos(a;n) cos{ujm) cos^{uJo)dr 

and similar terms with a different cosine factor replaced by a sine factor and corresponding 
prefactor. The product of the six trigonometric functions is odd about the point 7r/4. Thus 
this product integrates to 0 over each period. On the period [R + 2'Kk, R + 27r{k + 1)) with 
any nonnegative integer k, we may thus replace the weight r~^ by the difference between 
r~^ and its mean over that interval. This difference is bounded by 6r“^7r on that interval, 
hence we may estimate the sum of terms arising from the second main term of Corollary 

[Hby ^ 

37r((37)2 + (19)2 + (18)2+ 3)y 0^ r“^dr < 2.1 x 10“^^ . 

The sum of the six terms arising from the error terms of Corollary[^can be further estimated 
by ^ 

^((37)"^ + (19)^ +(18)^ +3) y 0^ < 1.64 X 10"^ . 

Next come fifteen terms of the original 2® — 1 error terms which have four main terms 
and two error terms of Corollary These benefit from an integration of the negative fourth 
power of r, and can be estimated by 

[372 X 192 + 372 X 182 + 192 X 18^ + 12 X 362] J ^ 2 ^ ^ ^ ^ ^ 

The remaining 2® — 1 — 6 —15 = 42 terms benefit from an integration of at least the negative 
fifth power of r, and are estimated even more crudely as 

42 X [372 X 192 X 182] J r-^dr < 4.5 x 10 "^° . 

Adding all these error contributions yields 


4((n + m)2 — 1) 


pOi 

Jr 


00 ^ 2 

nr 




^2,high 


-Ir 


1 , high I ^ 5.5 X 10 


-9 
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We next turn to the low integrals. We recall the Newton-Coates rule 

[' fix) dx = F{f) 


with 


Fif) = 140-'(41/(0) + 216/(1) + 27/(2) + 272/(3) + 27/(4) + 216/(5) + 41/(6)), 

which is valid for all real polynomials / up to degree 7. For any eight times continuously 
differentiable function / on [0, 6], we have that 


(73) 


fix) dx - Fif) 


6 ^ 

< — sup 
0 5e[o,i] 




8 ! 


A well-known argument shows that polynomials of degree eight extremize this inequality. 
It is then a straightforward matter of checking that polynomials of degree eight, whose 


eighth derivative is constant, realize the optimal constant ^ promised by (|73). 


Now let Fa^w be the suitably scaled and translated Newton-Coates formula which inte¬ 
grates polynomials of degree 7 on the interval [a, a + btc] exactly. Then, by rescaling, 


ca+6iti 


fix)dx - Fa,wif) 


6 ^ 

< W — 


sup 

5e[a,a+6ui] 


l/">K)l 


81 


Now assume that the length of the interval [a, b] is an integer multiple of brc, say 6wN. 
Then partitioning this interval into N intervals of length Qw and applying the Newton- 
Coates formula on each interval yields 


rb ^ ) 03 

/ fix)dx — E Fa+kw,wif)\ < ib-a)w^— sup 

Jo. k=0 


.63 |/( 8 )( 0 | 

?e[a,b] 


8 ! 


We cut the interval [0,7?] into [0,5"] U [5,7?] with 5 = 3600. On the interval [0,5], we 
estimate the eighth derivative of the functions 

fir) = Jn+mir)Jnir)Jmir)JQir)r 

and 

fir) = Jn+mir)Jnir)Jmir)Jiir)Joir)r 

using the Cauchy integral formula for the circle of radius 1 about r together with the trivial 


bound (10) to obtain the estimate 


|/(®)(r)| <8!e®(5 + l) . 

Approximating the integral over [0, 5] by the above summation rule with width w = 0.003 
gives the error bound 

63 

5(.003)®—e®(5 -b 1) < 1.49 x 10"^ . 













38 


DIOGO OLIVEIRA E SILVA AND CHRISTOPH THIELE 


On the interval [5, R], we estimate the eighth derivative of / again by the Cauchy integral 
formula with circles of radius one. We use the estimate from Corollary for to obtain, 
for >5 — 1 and 9‘(z) < 1, 

\Jn{^)\ ^ ( “Tt) (1 + cosh(9‘(z)) X 1.01 , 

\7r|2;| / 

and similarly for J~. Estimating the product of the various terms analogous to (l+n^/5) x 
1.01 by 3, we then obtain 

|/(®)(r)| < 3 X 8! X (cosh(l))6(i2+ 1) . 

Approximating the integral over [5, R] by the above summation rule with width w = 0.05 
gives the error bound 

3x{R- S)w^^(—J—-Y{cosh{l)f{R + 1) < 1.42 x 10“^ . 

5 V7r(5 — Ij/ 

Collecting error terms, we obtain 

I Jo - 1-2798 X 10"® - E[o, 5 ] - F[s,r] \ < 0.85 x 10"® 

if n is even, and 

I Jo - 0.256 X 10"® - F[o, 5 ] - F[s,r]\ < 0.85 x 10"® 

if n is odd, and 

\h - 0.256 X 10"® - F[o, 5 ] - F[s,r]\ < 0.85 x 10"® 
for any n, where and F^g^^ are the quadrature formulae described above for the 

corresponding integrals. 

We evaluate and F^g^^ using Mathematica. Products of Bessel functions at the 

grid points are computed with 20-digit precision, and the corresponding rounding errors 
for ^[ 0 ^ 5 ] can be safely estimated by 0.05 x 10“®. As an example, in the case n = 14 

and m = 4 for Jq, we use the following code to compute 


BJ[x_] := N[BesselJ[18, x] =t= BesselJ[14, x] * BesselJ[4, x] * BesselJ[0, x]® =t= x, 20] 

BJSA := 41 * BJ[0] + 

82 * Sum[BJ[x], {x, 18/1000, 3599982/1000,18/1000}] + 

216 * Sum[BJ[x], {x, 3/1000, 3599985/1000,18/1000}] + 

27 * Sum[BJ[x], {x, 6/1000, 3599988/1000,18/1000}] + 

272 * Sum[BJ[x], {x, 9/1000, 3599991/1000,18/1000}] -b 
27 * Sum[BJ[x], {x, 12/1000, 3599994/1000,18/1000}] -b 
216 * Sum[BJ[x], {x, 15/1000,3599997/1000,18/1000}] + 

41 * BJ[3600] 

(BJSA* .018)/840 





ESTIMATES FOR CERTAIN INTEGRALS OF PRODUCTS OF SIX BESSEL FUNCTIONS. 


39 


For m = 0 and even n, Table lists upper bounds for the quantities 
3 13 1 


n 3271^ (n — l)n(n + 1 ) 
on top of each entry, and for 
113 1 

r- + 


- 1.2798 X 10"® - F[o,5 ] - Fi 


[S,R] 


+ 0.9 X 10"^) 


lOOn^ 


- 0.256 X 10-6 - F[o,s] - 


+ 0.9 X 10““ lOOn 


dvr^ n 3271^ (n — l)n(n + 1 ) 

at the bottom of each entry, with the appropriate quadrature formulae ^[ 0 ^ 5 ] and 
described above. For m = 0 and odd n, it similarly lists upper bounds for 
3 13 1 


dvr^ n 3271^ (n — l)n(n + 1 ) 


- 0.256 X 10 


-6 


“ -^[0,5] “ ^[S,R] 


+ 0.9 X 10” 


lOOn^ 


on top, and for 


/ 1 1 3 1 

V dvr^ n 3271^ (n — l)n(n + 1) 


- 0.256 X 10-6 - F[o,5] - F^s,R] 


+ 0.9 X 10-®) 


lOOn"^ 


at the bottom. Thus each entry on the first column {m = 0) of Table divided by 100, 
provides a constant c for which the estimate of Theorem holds for the corresponding n 
with c in place of 0.002 or 0.0015. The entries of Table for m > 0 are analogous. 

The poorer constants near n = 19 are artificial and due to the chosen numerical accuracy 
0.9 X 10-®; note that, for this value of n, the quantity 0.9 x 10-®n^ is already close to O.OOld. 
The very good constants at m = d are due to the extra term in the expansion that has 
been elaborated in that case. 
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n\m 0 2 4 6 8 10 12 14 16 18 

2 ^^l4 
.64 .03 

3 ^^16 
.21 .05 

4 .33 .16 .03 

.16 .04 .01 

5 .26 .15 .02 

.12 .04 .01 

6 .22 .15 .02 .11 

.10 .04 .01 .06 

7 .19 .14 .02 .09 

.09 .04 .01 .05 

8 .17 .13 .02 .08 .02 

.08 .04 .01 .05 .02 

9 .15 .13 .02 .07 .02 

.07 .04 .01 .05 .02 

To .14 .13 .02 .07 .02 .02 

.07 .04 .02 .04 .02 .02 

11 .13 .13 .02 .07 .03 .02 

.07 .04 .02 .04 .02 .02 

12 .13 .13 .03 .07 .03 .03 .03 

.07 .05 .03 .05 .03 .03 .03 

13 .13 .13 .04 .07 .04 .03 .03 

.08 .05 .03 .05 .04 .03 .03 

14 .13 .13 .05 .07 .05 .04 .04 .04 

.08 .06 .04 .06 .05 .04 .04 .04 

15 .14 .14 .06 .08 .06 .06 .06 .06 

.09 .07 .06 .07 .06 .06 .06 .06 

16 ■ .15 .15 .07 .09 .07 .07 .07 .07 .07 

.10 .09 .07 .08 .07 .07 .07 .07 .07 

17 .16 .17 .09 .11 .09 .09 .09 .09 .09 

.12 .10 .09 .10 .09 .09 .09 .09 .09 

18 .18 .19 .11 .13 .11 .11 .11 .11 .11 .11 

.14 .13 .11 .12 .11 .11 .11 .11 .11 .11 

19 .20 .20 .14 .15 .14 .14 .14 .14 .14 .14 

.16 .15 .14 .14 .14 .14 .14 .14 .14 .14 

Table 1 
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